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ABSTRACT 

Cosmic shear has been identified as the method with the most potential to constrain 
dark energy. To capitahse on this potential it is necessary to measure galaxy shapes 
with great accuracy, which in turn requires a detailed model for the image blurring 
by the telescope and atmosphere, the Point Spread Function (PSF). In general the 
PSF varies with wavelength and therefore the PSF integrated over an observing filter 
depends on the spectrum of the object. For a typical galaxy the spectrum varies across 
the galaxy image, thus the PSF depends on the position within the image. We estimate 
the bias on the shear due to such colour gradients by modelling galaxies using two 
co-centered, co-elliptical Sersic profiles, each with a different spectrum. We estimate 
the effect of ignoring colour gradients and find the shear bias from a single galaxy 
can be very large depending on the properties of the galaxy. We find that halving the 
filter width reduces the shear bias by a factor of about 5. We show that, to first order, 
tomographic cosmic shear two point statistics depend on the mean shear bias over the 
galaxy population at a given redshift. For a single broad filter, and averaging over a 
small galaxy catalogue from Simard et al. (2002), we find a mean shear bias which 
is subdominant to the predicted statistical errors for future cosmic shear surveys. 
However, the true mean shear bias may exceed the statistical errors, depending on 
how accurately the catalogue represents the observed distribution of galaxies in the 
cosmic shear survey. We then investigate the bias on the shear for two-filter imaging 
and find that the bias is reduced by at least an order of magnitude. Lastly, we find 
that it is possible to calibrate galaxies for which colour gradients were ignored using 
two-filter imaging of a fair sample of noisy galaxies, if the galaxy model is known. For 
a signal-to-noise of 25 the number of galaxies required in each tomographic redshift 
bin is of order IC*. 

Key words: galaxy shapes 



1 INTRODUCTION 

Cosmic shear is the weak distortion of distant galaxy images 
due to the weak gravitational lensing of light by intervening 
matter. Light rays emitted by nearby galaxies follow similar 
paths through the Universe and thus their images are coher- 
ently distorted. A statistical analysis of correlations in the 



cosmic shear signal on different scales therefore provides cru- 
cial information about the distribution of dark matter and 
thus the cosmological model. Several upcoming and future 
observational surveys plan to capitalise on the potential of 
cosmic shear for discovering the nature of dark energy. These 
include the ground-based projects the Kilo-Degree Survey 
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(KIDS), Pan-STARRS B Suftara Measurement of Images 
and Redshifts (SuMIRe)Bthe Dark Energy Survey (DES)0 
and the Large Synoptic Survey Telescope (LSST)0, and the 
proposed space missions Eucli(^ and lyi^/i^Sllj. 

As surveys become more ambitious and the pro- 
jected statistical uncertainties on cosmological parameters 
reduce, the biases associated with shear measurement 
must be understood and controlled with increasing ac- 
curacy. Considerable work has been done on the main 



images (see Hevmans et al. 20061: 


Massev et al.l l2007l: 


Bridle et al.ll200g|.l2010l:lKitchinE et al. 


I2OI0I. and references 


therein), requirements on measurement and calibration 
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2005; 
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Kirk et al.ll2010l). 
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In this paper we consider a crucial aspect of shear 
measurement: deconvolution of the Point Spread Function 
(PSF). As light from a galaxy passes through the atmo- 
sphere and telescope optics and onto the detector it is con- 
volved with a kernel known as the PSF. The PSF is gener- 
ally not circular and we typically want to measure cosmic 
shear for galaxies down to the size of the PSF. The PSF 
must be accurately determined either from the galaxy spec- 
trum and a detailed model of the telescope and/or from stars 
in the image, which can be treated as point objects before 
the convolution. Several papers have put forward methods 
for quantifying the spatial (across the detector) and tem- 
poral variations in the shape of the PSF (e.g. Ijain et al] 
I2OO6I : lAmara et al.ll201ol ). as well as looking at the number 
of noisy stars needed to accurately ca librate the PSF model 
l|Paulin-Henriksson et aLll2008l . |2009| ) . 

A further issue is the wavelength dependence of the PSF 
which becomes more important when galaxies are imaged in 
wide filters. The PSF is a function of wavelength and there- 
fore the observed image consists of the galaxy image at each 
wavelength, convolved with the PSF at e ach wavelength, 
and th en integrated over the filter response. ICvpriano et al.l 
investigated the first order impact of galaxy spectral 
energy distributions (SEDs) combined with a wavelength 
dependent PSF. If the spectrum of the galaxy is not known 
then in principle this makes it impossible to deconvolve the 
PSF, even if the PSF shape as a function of wavelength 
is perfectly known. ICvpriano et al.l (|201(j ) have shown that 
the colour information required from ground-based surveys, 
such as the DBS, for estimating galaxy photometric redshifts 
is enough to provide the required spectral information, if 
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galaxies have the same spectrum at all points on their im- 
age. 

In this paper we investigate the second order effect of 
galaxy internal colour gradients. This has been hig hlighted 
in th e context of upcoming space missions (Cvpriano et al.l 
|2010| ). and Hirata & Bernstein (priv. comm.) have shown 
that, if ignored, the effect can produce shear measurement 
biases more than an order of magnitude larger than the re- 
quirements for upcoming surveys. 

Galaxies are often described by two-component models 
in which the bulge component is redder then the disk com- 
ponent. In this case there is no single PSF across the im- 
age if the PSF is a function of wavelength. Using one SED 
for the bulge component and another for the disk we simu- 
late galaxy images and recover shear values using the PSF 
for the composite galajcy spectrum across the whole image. 
The resulting bias is compared to the requirement for future 
cosmic shear surveys to determine whether a single optical 
imaging filter is sufficient if galaxy color gradients are ig- 
nored. We use a catalogue of galaxies to find a reasonable 
range of galaxy parameters. We then repeat the calculation 
with two filters, allowing estimates of the bulge and disk 
spectra separately. Finally, we investigate the feasibility of 
calibrating the bias on a subset of galaxies imaged in more 
than one filter. 

The paper is organised as follows. In Section [2] we de- 
scribe the galaxy and telescope models used in the simu- 
lations. In Section [3] we briefly outline the equations de- 
scribing the distortion to galaxy images from weak gravita- 
tional shear. We then describe in detail the simulations and 
method used to estimate the biases on the shear and set 
out the requirements on the bias for future surveys. In Sec- 
tion[4]we quantify the bias on shear estimation from internal 
colour gradients for single filter imaging. In Section [5] we in- 
vestigate the improvement gained by imaging in two filters. 
We then quantify the number of galaxies needed to calibrate 
the biases in Section [6] Finally we discuss the implications 
of these results in Section [T] 



2 THE GALAXY AND TELESCOPE MODELS 

Here we describe the models used to generate the galaxy 
and PSF images used in the simulations. We also outline 
the procedure for convolving the sheared multi-component 
galaxy images with the PSF. 



2.1 The telescope model 

The image of a point source is broadened as a result of 
two dominant wavelength-dependent processes: (1) diffrac- 
tion giving rise to an Airy disk with size proportional to 
wavelength and (2) the detector modulation transfer func- 
tion (MTF) which tends to spread out higher energy photons 
more than lower energy photons. The PSF intensity map as 
a function of x for a galaxy with SED S{X) is given by 



J^°° IAx-X)S{\)T{\)dX 
J^°° S{\)T{X)dX ' 



(1) 



where Jp(a;;A) is the normalised PSF map at each wave- 
length, T(A) is the instrumental plus filter response and 
S{X)T{X)dX is the total observed flux from the galaxy. 
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Figure 1. Disk (upper blue solid curves) and bulge (lower red 
solid curves) for a CWW-Sbc (left panels) and CWW-Im (right 
panels) disk spectrum for a galaxy at redshift z = 0.6 (upper pan- 
els) and z = 1.2 (lower panels). The bulge spectrum is CWW-E. 
The response functions T(A) are shown for the F\ (black dotted) 
and F2 (green dash-dotted) filters. The dashed lines show the 
linear fits to the fluxes in the -Fi and filters. 



Throughout the paper we assume a top-hat response 
function with r(A) = 1. For simplicity we model the PSF 
intensity at each wavelength, /p(x;A), by a Gaussian with 
'half-light radius' (radius enclosing half the total flux) given 
by 

rp(A) = rp.o (^j , (2) 

where Ao = 520nm and rp,o = 0.7 pixels. The PSF param- 
eters Ao and rp,o are chosen so that the Full Width at Half 
Maximum (FWHM) of the flux-weighted PSF of the fidu- 
cial galaxy (see Section[2]4]) is approximately 1.7 pixels. The 
PSF ellipticity and orientation are fixed throughout the pa- 
per at ep = 0.05 and <^p = 0° respectively. The FWHM and 
ellipticity of the PSF are chosen to approximately represent 
the values expected in future cosmic shear surveys. 

We have used an elliptical Gaussian throughout the pa- 
per to model the PSF at each wavelength. For more realistic 
estimates a PSF similar to the instrument in question should 
be used, for example an Airy disk for a space mission. 

The fiducial filters we consider in this paper are a broad 
top-hat filter with width 550-900 nm, and a narrower filter 
with width 725-900 nm (hereafter F\ and F-i respectively). 
We choose F\ to correspond to the filter currently under 
consideration for the Euclid satellite. F2 is half the linear 
width of the Fi filter and at the red end of the Fi filter. 
The filter response functions are shown in Fig. [1] along with 
example bulge and disk galaxy spectra, which are introduced 
in Section [231 



2.2 The galaxy model 

In this paper we simulate two-component galaxies with a 
realistic range of profile shapes and colour gradients. The 
bulge and disk are modelled by Sersic profiles (Sersic 1968). 
The Sersic intensity at position x — (x, y) is given by 



where xq is the centre, Iq is the peak intensity, Us is the 
Sersic index and C (proportional to the inverse covariance 
matrix if ris = 0.5) has elements 
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where is the angle (measured anti-clockwise) between the 
X-axis and the major axis of the ellipse and the minor to 
major axis ratio is b/a. 

The Sersic index defines the profile 'type', with ris = 
0.5, 1 and 4 for Gaussian, exponential and de Vaucouleurs 
profiles respectively. If k is defined as fc = 1.9992ns — 0.3271 
then for a circular profile ro = a = 6 is the half-light ra- 
dius (Note that for a Gaussian profile and are the 
2D variances if fc = 0.5; for the circular exponential profile 
h — a = b is known as the 'scale length' when fe = 1.) 

The fiducial model used in this paper for the bulge is 
a de Vaucouleurs profile and for the disk an exponential 
profile. The FWHM is related to the half-light radius (for a 
circular profile) via 



FWHM = 2r-< 



(7) 



The total flux (integrated to inflnity) emitted by a galaxy 
described by a Sersic proflle with index ris is given by 



F = 2nnsk~ 



^r,'r(2ns)/o 



(8) 



where F is the gamma function. 



Ig(x) = Joe 



-k[(x-Xo)^ Cix-Xo)]" 



(3) 



2.3 The galaxy catalogue 

The results will depend on the exact assumptions made 
about the galaxy properties and therefore we use a galaxy 
catalogue to estimate the average qu antities required. The 
catalogue we use is from a study by ISimard et all (|2002[ ) 
of Hubble Space Telescope (HST) WFPC2 F606W (470- 
720 nm; 'F-band') and F814W (708-959 nm; '/-band') ob- 
servations of the 'Groth Strip'. The catalogue lists galaxy 
photometric structural parameters obtained from a two- 
dimensional bulge plus disk decomposition of the galaxy sur- 
face brightness profiles. The bulge profile is modelled as a 
de Vaucouleurs and the disk as an exponential. 

ISimard et al.l l|2002l ) produce two sets of catalogues, one 
in which fits are performed separately in each filter, and 
a second catalogue in which the fits are performed simul- 
taneously in each bandpass. We use the second catalogue, 
which is from the DEEPl (Deep Extragalactic Evolutionary 
Probe) Keck LRIS spec troscopic survey of the Groth Strip 
l|Vogt N. P. et al.|[2005l ). The sample is magnitude-limited 
and contai ns several objects w hich were prioritized for se- 
lection (see lWeiner et al]|2005l ). The catalogue contains 632 
galaxies in the magnitude (Vega) range {V + 1)/2 < 24 and 
the median redshift is 0.65. 

The catalogue is well-suited to our study because the 
broad V and / bands used in the HST observations cover 



© 2008 RAS, MNRAS 000.[llfT6l 



4 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

redshift redshift redshift 

Figure 2. Bulge colour (left), disk colour (middle) and bulge minus disk colour (right) as a function of redsiiift for a CWW-Sbc (blue 
solid) and CWW-Im (green dashed) disk spectrum. The bulge spectrum is CWW-E (red solid). Fluxes are measured in the V and / bands. 
The shaded regions show the central 68 per cent of the colour distribution measured from galaxies in the catalogue (see Section 12.31 1 . 



a similar wavelength range to the broad Fi filter. The cata- 
logue provides bulge and disk scale radii (fixed at the same 
value in each filter) and the magnitudes of the bulge and the 
disk in each filter. 

We compute linear approximations to the bulge and 
disk spectra, S{X), across the V and / bands from the fluxes 
in each filter using the AB magnitudes given in the cata- 
logu4ll- The bulge SED between the midpoints of the V and 
/-bands (595-833.5 nm) is given by ^(A) — rribA + Cb, where 
rrib and Cb are the gradient and y-intercept of the spectrum 
respectively, and similarly for the disk. The spectrum is as- 
sumed to be flat between 470-595 nm and 833.5-959 nm, 
such that ^(A) = m]jXv^i^ + "^t in the first half of the V- 
band and S{X) — m\^Xi^.^ + Cb in the second half of the 
/-band, where Xv^^^ = 595nm and Xi^.^ — 833. 5nm. The 
gradient and intercept of the bulge and disk spectra are com- 
puted so that the total bulge and disk fluxes in the V and / 
bands are equal to the values in the catalogue. If S{X) goes 
negative in either band then the spectrum is re-set to be flat 
with a different value on either side of 720 nm, with the con- 
straint that the flux between 470-720 nm is equal to the flux 
in the l/-band, and the flux between 720-959 nm is equal to 
the flux in the /-band minus the flux between 708-720 nm 
calculated from the flat spectrum between 470-720 nm. 

The exact distribution of galaxies in a particular sur- 
vey will depend on the parameters of the survey (e.g. the 
median redshift), and this must be taken into account when 
using the results obtained from the catalogue described here 
to make predictions for future cosmic shear surveys. We in- 
vestigate the sensitivity of the results to the galax;y redshift 
and colour distributions in the catalogue in Section [4.21 

We also note that the photometric and structural pa- 
rameters obtained in the ISimard et al.l (i2002) study will 
not precisely represent the underlying distribution of galaxy 
colours and shapes in the DEEPl survey. This is because (i) 
the study fits a simple two-component model to the galaxy 
images, when in fact galaxy morphologies may be complex 
and (ii) galaxy images are noisy wh ich may bias the fitted 
parameters (see lHaussler et al]|2007l ). 



m = — 2.51ogio(F) — 48.6, where m is the AB magnitude and 
F is the flux. 



2.4 The fiducial galaxy model 

We also consider a fiducial galaxy model which is the 'aver- 
age' galaxy from the catalogue. The fiducial galaxy param- 
eters are ns,b = 4, ris.d = 1, rc,b/rc,d = 1-1, B/T = 0.25 
and z = 0.9. The bulge and disk ellipticities are set equal 
to eh = ed = eg — 0.2. The fiducial values chosen for the 
ratio of the bulge-to-disk half light radii and the bulge-to- 
total flux ratio are the mean values from the catalogue. The 
central 68% of the rc,h/rc,d distribution is approximately 
l.llg g. The flducial redshift (z — 0.9) is the median value 
for Euclid. 

We use a CWW-E spectrum for the bulge throughout 
the paper, and e ither a CWW-Sbc or CWW-Im spectrum 
for the disk (see lColeman et al.|[l980h . We use the CWW- 
Sbc spectrum for the fiducial disk spectrum. The spectra 
observed in the wavelength range 500-1000 nm are shown 
in Fig. [T] for two different source redshifts. We show in Fig. [5] 
that these spectra provide a good representation of the dis- 
tribution of bulge and disk colours found in the catalogue, 
with the CWW-Sbc spectrum representing the disk V — / 
colours more closely than the CWW-Im spectrum. 

The bulge and disk half-light radii are rescaled (keeping 
the ratio of the bulge half-light radius to disk half-light ra- 
dius constant at the fiducial value) throughout the paper so 
that the ratio of the PSF-convolved galaxy FWHM to PSF 
FWHM, rgp, is fixed at 1.4 for a circular PSF and galaxy. 



2.5 Convolution 

The PSF-convolved galaxy intensity at x is given by 

/•oo 

/g,p(x)= / /g(x;A)*/p(x;A)r(A)dA (9) 
Jo 

where /g(a;;A) is the galaxy intensity per unit wavelength. 
For multiple component galaxies with components i this can 
be written as 



i 

X,°°Si(A)T(A)/p(a;; A)dA 



/p,i(a;) = 



(10) 
(11) 



/;°Si(A)T(A)dA 
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where /g,i(a;) is the galaxy intensity integrated over all wave- 
lengths and Si{X)T{X)d\ is the total observed flux from the 
ith component per unit wavelength. 



3 QUANTIFYING THE BIAS ON SHEAR 
MEASUREMENTS 

We quantify the bias on cosmic shear measurements 
from using an incorrect PSF model for a given set of 
ga laxy parameters. We us e the 'r ing-test' m ethod described 



iNakaiima fc BernsteinI (|2007l ') (see also IVoigt fc Bridle! 
|2010D . We first briefly outline the equations describing the 
distortions to galaxy images from gravitational shear before 
describing the simulations in more detail. 



3.1 Gravitational Shear 

The two-component gravitational shear of a galaxy, given 
by 



7 = 71 + *72, 



(12) 



is related to the projected dimensionless gravitational po- 
tential ?/) of a thin lens between the source and the observer 
such that 
1 



71 



{ipll — i>22) , 72 = ^12 = V'21 



(13) 



where ipij — d^ijj /dOidO j and is the observed position of 
the source (see, for e.g.. iBartelmann fc Schneider 2001). In 
the weak lensing regime the reduced shear, g, is approxi- 
mately equal to the shear 7 where 



5 = 7/(1 - 1^), 



(14) 



and the convergence k is proportional to the projected sur- 
face meiss density of the lens. A pre-shear elliptical isophote 
galaxy image with minor to major axis ratio b/a and ori- 
entation of the major axis anticlockwise from the positive 
a;-axis, (j>, has a complex ellipticity 

e= = ^e^^^^ (15) 
a + b 

and is sheared to an ellipse with complex observed ellipticity 
e° given by 

e° = -£l±^ (16) 

l|Seitz fc Schneider! [l997h if K < 1. We will assume k << 1 
throughout this paper. 

3.2 The Simulations 

We simulate two-component galaxies at different orienta- 
tions with co-centered, co-elliptical profiles. Each galaxy is 
sheared by the same amount g. The estimated shear for each 
galaxy is the measured observed ellipticity. We use a ring- 
test to quantify the bias on the estimated shear for a given 
galaxy model and PSF estimate in the limit of an infinite 
number of orientations of the pre-sheared galaxy. For a per- 
fect shear measurement method the average observed ellip- 
ticity is then equal to the true input shear. In practice we 
can reduce the number of galaxies needed in the ring-test by 
simulating pairs of galaxies separated by 90 degrees (as in 



iMassevllioOTh . It can be shown from Eqn.[T6]that even for a 
perfect shear measurement method we need at least 3 pairs 
of galaxies in the ring-test to reduce the bias on the esti- 
mated shear to a negligible level. If the shear measurement 
method is imperfect (e.g. because the PSF model is incor- 
rect) then we may need more than 3 pairs to obtain a shear 
value that has converged to the value we would measure for 
an infinite number of galaxy orientations. We find, however, 
that 3 pairs of galaxies is sufficient to reach convergence. We 
check that doubling the number of pairs of galaxies does not 
change the results. 

PSF-convolved galaxy images are produced on postage 
stamps IS'^ pixels in size using the following procedure. The 
PSF and galaxy models are described in Sections l2.1l and l2.2l 
respectively. The complex ellipticities of the post-shear bulge 
and disk are computed using Eqn.[T6] The size of each com- 
ponent (r^ b S'lid fo,d) is rescaled by a factor (1 — |5p)~^ after 
shearing. The PSF convolutions are performed on a high res- 
olution image which is 17^ observational pixels in area and 
has each observational pixel divided into 7^ subpixels. We 
find that increasing the number of subpixels used in the con- 
volution beyond 7^ has a negligible impact on the results. 
We find that it is necessary to create the galaxy image at an 
even higher resolution before binning up to make the galaxy 
image used in the convolution. We find it is sufficient to do 
this only on the central 5^ subpixels of the convolution im- 
age, but that each of these subpixels must be divided into 
25^ subsubpixels. The intensity in each subpixel is found by 
calculating the intensity at the centre of each subsubpixel 
and summing all subsubpixels within each pixel. The inten- 
sity in the remaining subpixels of the image is assumed to 
be equal to the value at the centre of that subpixel. Integra- 
tion of the intensity within each of the central 5^ subpixels 
of the image is necessary for highly peaked profiles, such as 
the de Vaucouleurs profile used to represent the bulge sur- 
face brightness distribution. The convolution with the PSF 
is performed using FFTs using sufficient padding such that 
the result of the convolution is identical to that from per- 
forming the convolution in real space. The high resolution 
grid is larger than the postage stamp to allow for light scat- 
tered back into the postage stamp due to convolution with 
the PSF. We check that increasing the resolution parameters 
does not change the results. 

The high resolution PSF-convolved bulge and disk im- 
ages are then binned up by a factor of seven and cut down 
to produce an image on a postage stamp 15^ pixels in size. 
The observed galaxy image is the linear sum of the bulge 
and disk PSF-convolved images. The bulge and disk are co- 
centric and co-elliptical and the peak of the intensity profile 
is within the central pixel of the image grid. The fiducial 
centre used is {xo,yo} = {0.1,0.3} pixels relative to the 
centre of the postage stamp for the galaxy aligned along the 
X-axis. (In the ring-test the galaxy is rotated around the 
centre of the grid). The PSF FWHM is 1.7 pixels and the 
FWHM of the PSF-convolved galaxy image is 1.4 times the 
PSF FWHM. 

The best-fit parameters of the bulge and disk compo- 
nents are computed on the image (low resolution) grid for 
each angle in the ring test by minimising the between the 
true PSF-convolved galaxy image and the PSF-convolved 
galaxy image using the estimated PSF model. This process 
is repeated for two different input shears: {71,72} = {0,0} 
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and {7!;, 72} ~ {0.01,0.02}, where the superscript 't' refers 
to the true value of the shear. For each input shear we obtain 
a shear estimate, 7. We quantify the bias on the shear esti- 
mator in terms of multip hcative and add itive errors, rui and 
Ci respectively, following iHevmans et al^ I iooei), such that 



7i = (1 + mi)7- + d 



(17) 



where the subscript i refers to the two shear components 
and we assume there is no cross contamination of e.g. 71 
depending on the value of 72 or vica versa. 



3.3 Survey requirements 

We use the requirements on the multiplicative and additive 
sh ear biases for cosmi c shea r two poi nt statistics compute d 
in lAmara k Refregied l|2008l ) (see also lKitching et al.ll2009l ). 
These are set so that the systematic error is equal to the 
statistical error, where the statistical uncertainties depend 
on the survey area, depth and galaxy number density. For 
a medium-deep survey of 20,000 square degrees (hereafter 
^Euclid-like survey') the upper limits on the total contribu- 
tion to the multiplicative and additive biases are 10~^ and 
3 X 10"'* respectively. Biases a factor of 1-2, 2-5 and 5- 
20 below the Euclid-like survey requirements are shown in 
the figures by light, medium and dark grey shaded regions 
respectively. 

The requirements on the she ar bias for a particular cos - 
mic shear survey as computed in I Amara fc Refregiej (|2008| ) 
must encompass all systematics involved in the survey (e.g. 
shear measurement, charge transfer inefficiency or CTE, in- 
trinsic alignments etc.), and thus to be safe for a Euclid-like 
survey we require the bias from the colour dependence of 
the PSF to be of order a factor of 10 below the upper limits 
on the multiplicative and additive biases stated above. 



4 SINGLE FILTER IMAGING 

We calculate shear measurement biases for imaging in a sin- 
gle filter, first considering the dependence of the bias on 
filter width for some example galaxies, and then calculating 
the mean bias for single broad-filter imaging by averaging 
over the catalogue. 



4.1 Dependence on filter width and galaxy 
properties 

We first compute the multiplicative and additive biases on 
the shear for a two-component galaxy imaged in a single visi- 
ble filter of increasing width with central wavelength 725nm. 
For single filter imaging we have no resolved colour informa- 
tion; we assume however that the galaxy spectrum from the 
composite bulge plus disk image is known perfectly from un- 
resolved observations in multiple filters. Such observations 
are already required for estimating galaxy photometric red- 
shifts. The PSF model is described in Section 12.11 We use 
the PSF model obtained using the SED from the composite 
galaxy spectrum (equal to the flux- weighted sum of the true 
bulge and disk PSF models) for both the bulge and disk 
PSFs. The galaxy model is described in Section 12.21 The 
bulge-to-total flux ratio is kept constant as the filter width 
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Figure 3. Absolute value of the multiplicative (top) and addi- 
tive (bottom) shear bias as a function of filter width for imaging 
in a single filter. The filter transmission function is a top-hat 
and the central wavelength of the filter is 725 nm. Open squares 
(crosses) show mi,ci (m2,C2). The fiducial galaxy parameters are 
2 = 0.9,re,bAe,d = 1.1, -B/T = 0.25, a CWW-Sbc disk spectrum 
and a CWW-E bulge spectrum (red solid). Results are also shown 
for 2 = 1.4 (blue dashed), b/rj, ^ = 0.4 (green dot-dashed) and 
a CWW-Im disk spectrum (black dotted), with all other parame- 
ters kept at the fiducial values above. Note that C2 is zero because 
the PSF is aligned along the x-axis. The shaded regions from light 
to dark indicate where rrii is a factor of 1-2, 2-5 and 5-20 below 
the Euclid requirement on the bias. 



increases. Variations in the shear bias values are thus dom- 
inated by the increase in filter width, and not by a change 
in the galaxy shape. 

The biases are computed using the ring-test method 
described in Section 13.21 We assume the galaxy model is 
known. There are seven free parameters in the fits to the 
single filter galaxy images: xo,yo,eg,(j)g,re,h,re,d and 7o,b, 
where /o,b is the peak intensity of the bulge profile. Note 
that the peak intensity of the disk profile is known from 
'"e.b, re,d, lQ,h and the total galaxy flux. We assume the total 
galaxy flux would be obtained from other observations such 
as the unresolved observations used to obtain the composite 
spectrum. 

Results are shown in Fig. |3] for the fiducial galaxy (red 
solid lines), as well as for a CWW-Im disk spectrum (black 
dotted lines) and for a different value of the galaxy redshift 
(blue dashed lines) and r^b/rcd ratio (green dot-dashed 
lines) . The shear biases decrease as the filter width decreases 
because the bulge and disk spectra are in general more sim- 
ilar for smaller widths (and identical for vanishingly small 
widths). For surveys using wide-band imaging with a filter 
width of 350 nm (the Fi filter) the multiplicative biases for 
the fiducial galaxy considered here are a factor of 3 below the 
EucUd-like survey requirement. The additive bias is a factor 
of 20 below the requirement. The size of the bias is strongly 
affected by the ratio of the bulge to disk half-light radius, 
as well as the internal colour gradient, with the multiplica- 
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Figure 4. Absolute value of the multiplicative shear bias for 
imaging in the _Fi filter as a function of the galaxy parameter 
(clockwise from top left): bulge Sersic index, bulge-to-total flux 
ratio, galaxy ellipticity, position of peak intensity relative to the 
centre of the postage stamp (for the galaxy aligned along the 
X-axis; see Section 13.21 xq = 0.1). All other galaxy parameters 
are fixed at the fiducial values (see Section 12.41 1. The true bulge 
and disk spectra are CWW-E and CWW-Sbc respectively. The 
estimated PSF model for the bulge and the disk components is the 
model from the composite galaxy spectrum. Solid open squares 
(dashed crosses) show mi (m2 ) . Shaded regions as in Fig. |3] 

tive bias increasing to a factor of ~4 above the requirement 
and additive bias a factor of ~2 below the requirement for 
the single departure from fiducial of Tch/rcd ~ 0.4. We 
note that the distribution of bulge to disk half-light radii 
in the galaxy catalogue described in Section 12.31 has mean 
?'e,b/''G,d = 1.1, and the central 68% of the distribution is in 
the range r-o,b/f'e,d = 0.3 — 1.6. 

If we simultaneously make the above multiple depar- 
tures from the fiducial galaxy the biases are compounded, 
resulting in a multiplicative bias 30-50 times larger than 
shown by the red solid curve, depending on the filter width 
(e.g. for a filter width of 350 nm, m ~ 10~^ for z — 1.4, 
rc,h/rc,d = 0.4, a CWW-Im disk spectrum and a CWW-E 
bulge spectrum). For our co-elliptical simulations the mul- 
tiplicative biases dominate the additive biases and thus we 
consider only the multiplicative biases in the rest of the pa- 
per. We note that for the fiducial galaxy parameters the 
multiplicative bias increases (decreases) at each filter width 
by a factor of about 1.5 when the PSF ellipticity is doubled 
(halved) from the fiducial value of Cp = 0.05. The additive 
bias increases (decreases) at each filter width by a factor of 
about 3 when the PSF ellipticity is doubled (halved) from 
the fiducial value. 

We now show how different galaxy parameters affect 
the bias. In Fig. |4] we plot the multiplicative shear biases 
Eis a function of bulge Sersic index, bulge-to-total flux ratio, 
galaxy ellipticity and position of peak intensity within the 
central pixel of the postage stamp. We use the same fiducial 
values as in Fig. |3]with a filter width of 350nm (_Fi filter). 
The disk profile is an exponential. We see that changing 
these parameters within the ranges shown does not increase 
the bias above the Euclid-\ike survey requirement. The bias 
is little affected by the position of the galaxy. However, we 



see that to determine whether single broad filter imaging will 
give biased shear measurements we need to know the true 
distribution of galaxy shapes and colours in the Universe. 

4.2 Averages over the galaxy catalogue 

As shown in Section 14.11 the biases on the shear depend 
on the bulge and disk spectra, the galaxy redshift, the ra- 
tio of the bulge-to-disk half-light radii, the bulge-to-total 
flux ratio, the galaxy ellipticity and the bulge Sersic index. 
In this section we compute the biases for single broad-fllter 
imaging with the Fi fllter for a realistic range of galaxy pa- 
rameters obtained from a catalogue, described in Section [T3l 
The catalogue represents an observed range of galaxy colour 
gradients and bulge-to-disk half-light radii. The bulge Sersic 
index is Us.b = 4 and we fix the galaxy ellipticity at eg = 0.2 
for all galaxies in the catalogue. We see from Fig. |4] that by 
fixing the bulge Sersic index and galaxy ellipticity we may 
over- or under-estimate the mean bias by a factor of at most 
3, depending on the true distribution of these parameters in 
the Universe. 

For each galaxy in the catalogue we simulate separate 
'true' bulge and disk PSF images in the Fi filter by sub- 
stituting the linear spectra obtained from the bulge and 
disk magnitudes (see Section 12. 3p into Eqn [l] Bulge and 
disk galaxy images are computed and convolved with the 
PSF images as described in Section [3.21 The bulge and disk 
shape parameters re,b/f'e,d and B/T are taken from the cat- 
alogue. The bulge and disk half-light radii are scaled so that 
rgp = 1.4. In the fits to the 'true' PSF-convolved galaxy 
images we use the PSF image from the composite galaxy 
spectrum to estimate both the bulge and disk PSF images. 
We measure the bias on the shear using a ring-test (as de- 
scribed in Section [3. 2|) for each galaxy in the catalogue. We 
assume the galaxy model is known and fit for 7 free param- 
eters, described in Section [4. II 

In Fig. [5] we plot the multiplicative biases as a function 
of the galaxy internal V — I colour gradient (i.e. bulge V — I 
colour minus disk V — I colour). The bias scales approxi- 
mately linearly with the colour gradient; the shear is over- 
estimated for negative colour gradients (blue-cored galaxies) 
and under-estimated for positive colour gradients (red-cored 
galaxies). The mean bias (mi) over the 632 galaxies in the 
catalogue is ~ 3 x 10~*, a factor of about 3 below a Eu- 
clid-like survey requirement, and similar to the value ob- 
tained from the fiducial galaxy which was chosen to have 
average parameters from the catalogue. The bias does not 
change monotonically with colour gradient because the ratio 
of bulge to disk size is also changing within the catalogue. 

We now investigate the sensitivity of the mean bias 
to the colour distribution of galaxies in the catalogue. In 
Fig. [6] we plot the mean bias as a function of the percent- 
age of galaxies removed from the catalogue. We investigate 
the sensitivity to the colour gradient distribution by pro- 
gressively removing galaxies with the most negative colour 
gradient, and similarly for the most positive colour gradient 
(red dashed). The mean bias is equal to the mean bias over 
the whole catalogue for 0% removal. If we remove galaxies 
with the most negative colour gradients (blue-cored galaxies, 
shown by the blue solid line) then the mean bias increases 
relative to the mean bias over the complete catalogue by 
a factor of at most 3. If we progressively remove galaxies 
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Figure 5. Mean multiplicative bias (mi open squares; m2 
crosses) as a function oiV — I colour gradient using galaxies from 
the catalogue (see Section [2.3l l. Each colour gradient bin contains 
about 90 galaxies. The galaxy images are simulated in a single 
broad filter {Fi) with a top- hat transmission function. The esti- 
mated PSF model for the bulge and the disk components is the 
model from the composite galaxy spectrum. The error bars show 
the error on the mean bias in each bin (i.e. o-^^/y^n^, where 
is the standard deviation of the rrii distribution and rig the num- 
ber of galaxies in the bin). Vertical (horizontal) black dashed lines 
show where the colour gradient ((mi)) is zero. Shaded regions as 
in Fig. [3] 



with positive colour gradients (red-cored galaxies shown by 
the red dashed line), then the mean bias initially decreases 
before increasing to the Euclid-like requirement for removal 
of 60% or more of the catalogue. We note that removing 
galaxies with the most negative (or most positive) colour 
gradients increases the mean bias because the distribution 
of colour gradients in the catalogue is relatively symmet- 
ric about zero colour gradient, and the bias scales approxi- 
mately linearly with colour gradient. 

We also consider the sensitivity of the catalogue to the 
galaxy redshift distribution. For a Euclid-like survey with 
median redshift z = 0.9 the 4000A Balmer break is within 
the Fi filter, thus we expect to observe a significant fraction 
of red bulges. We plot the mean bias as a function of the per- 
centage of the 'blue-est' (most negative V — I colour) bulges 
removed from the catalogue. The mean bias increases to a 
factor of ~2 above the Euclid-like requirement for removal 
of 60% of the blue-est galaxies from the catalogue. Thus if 
red bulges are under-represented by the catalogue relative to 
the fraction expected in a Euclid-like survey the mean bias 
for single broad-filter imaging may be close to the statistical 
equals systematic error requirement of 1 x 10""^. 

As discussed in Section FS.Sl the mean bias needs to be a 
factor of 10 or more below the statistical equals systematic 
error requirement to be safe for a Euclid-like survey, and 
thus we consider options for reducing the bias below that 
for single broad-filter imaging. As shown in Section [4. II the 
mean bias from single broad-filter imaging may be reduced 
by imaging in a narrower filter. Alternatively, it may be 
possible to reduce the mean bias from single broad-filter 
imaging to the required level by imaging in an additional 
filter, or by calibrating the bias on a sub-set of galaxies. 




20 40 60 80 

Percentage of galaxies removed 

Figure 6. Absolute value of the mean multiplicative bias (mi 
open squares; m2 crosses) as a function of the percentage of galax- 
ies removed from the catalogue, starting with galaxies with the 
most negative colour gradient (blue solid), most positive colour 
gradient (red dashed) and most negative bulge colour (green dot- 
dashed). Shaded regions as in Fig. [3] 



5 TWO FILTER IMAGING 

In this section we investigate how far the biases on the shear 
may be reduced using two filter imaging. We repeat the ring- 
test calculation in Section 14.11 using a narrower filter, F2 
(725-900nm), in addition to the broad Fi filter. The nar- 
rower filter covers the second half of the wavelength range 
of the wider filter (see Fig. [T| . 

The 'true' PSF-convolved galaxy images are generated 
as described in Section|4Tl]using a CWW-E spectrum for the 
bulge and a CWW-Sbc or CWW-Im for the disk. For two 
filter imaging we estimate separate spectra for the bulge and 
disk components by fitting to the observed two-filter image 
and using information from the composite galaxy spectrum. 
We use different PSF models for each of the bulge and disk 
components (/p'^b ''■rid I^'^^ respectively), according to their 
estimated spectra. 

We first create linear spectra in which the bulge spec- 
trum is assumed to be flat between Afi,o and Afi,o.25 
and also from Afi.o.ts and Afi,i, where Xpij is the wave- 
length a linear fraction / across the Fi filter, for example 
Afi,o.25 = 637.5 nm. Therefore Sb{X) — mbAFi,o.25 + Cb in 
the first quarter of the Fi filter and ^^(A) = mbAFi.0.75 + Cb 
in the last quarter of the Fi filter (i.e. the second half of ^2), 
and similarly for the disk. The gradient and intercept of the 
bulge and disk spectra are re-computed at each iteration in 
the fit so that the estimated bulge and disk fluxes in Fi and 
F2 are equal to the integrated fluxes from the linear spec- 
tra. If Si{X) from the linear fits to the fluxes goes negative 
in either band then the spectrum is re-set to be flat with a 
different value on either side of the midpoint of Fi . Example 
linear fits to the bulge and disk fluxes are shown in Fig. [T]for 
redshifts 0.6 and 1.2. In both cases the bulge-to-total flux 
ratio is 0.25. 

Using just the linear spectra fitted to the two filter re- 
solved images is pessimistic because in reality we would have 
additional extra unresolved colour information used to pro- 
duce the composite spectrum. We therefore calculate the 
residuals between the PSF image obtained using the com- 
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Figure 7. Absolute value of the multiplicative shear bias, \mi\, as a function of the ratio of the bulge-to-disk half-light radii for CWW- 
Sbc (upper plots) and CWW-Im (lower plots) disk spectra for a galaxy at redshift 0.4 (blue solid squares), 0.8 (red dashed circles), 1.2 
(green dot-dashed crosses) and 1.6 (black dotted stars). Results are shown for imaging in one broad filter (_Fi; left-hand plots) and in 
two filters {Fi, F2; middle plots). The right-hand plots show the magnitude of the difference between the multiplicative bias obtained 
by imaging in one filter (-Fi) and the estimated bias from imaging in this filter using the bulge and disk shape parameters and spectra 
obtained from the two filter fits. The bulge-to-total flux ratio is B/T = 0.25 and the galaxy ellipticity is eg = 0.2. The bias values for 
two-filter imaging and a CWW-Im spectrum (lower middle plot) lie below the lower limit on the y-axis for z = 0.8. Shaded regions as in 
Fig. [3 



posite spectrum (equal to the flux-weighted sum of the true 
bulge and disk PSF images) and the flux-weighted sum of 
the bulge and disk PSF images obtained from the linear 
approximations to the spectra (/p_b and 7p respectively), 
equal to /p,res. The residuals are added to the bulge and disk 
PSF images obtained from the linear approximations to the 
spectra at each iteration in the flt such that the bulge PSF 
becomes Jp^J; = Jp'J, + /p,rcs, and similarly for the disk. 

The bulge PSF is then convolved with the bulge image 
and the disk PSF with the disk image and the best-fit galaxy 
parameters found by minimising the between the sum 
of these images and the true PSF-convolved galaxy image. 
There are eight free parameters in the fits, compared with 
seven free parameters in the single filter fits. The additional 
free parameter is the bulge amplitude in F2. (In practice we 
fit for the bulge amplitude in F2 and in the first half of -Fi 
to ensure that the fiux in _Fi is greater than the flux in F2). 
As with the single filter fits, the disk amplitude in each filter 
is determined from the bulge fiux in the filter and the total 
flux (assumed known). 

The multiplicative biases on the shear for the two fil- 
ter configuration are shown as a function of the ratio of the 
bulge-to-disk half-light radii for different galaxy redshifts in 
the middle panel of Fig. [T] For comparison the biases for the 
single broad-filter imaging are shown in the left-hand panel. 
The right-hand panel is described in Section [6] The central 



68% of the rcb/^cd distribution from the catalogue ranges 
from approximately 0.3 to 1.6. The figure shows that the ad- 
dition of a second narrower filter reduces the biases by more 
than an order of magnitude. We see that the biases for the 
two filter imaging are on average well within the Euclid-\ike 
survey requirements and thus two filter imaging on board 
a Euclid-like telescope would be sufficient to correct for the 
colour dependence of the PSF. 



6 BIAS CALIBRATION 

In this section we investigate the feasibility of calibrating 
the biases from internal colour gradients by imaging a sub- 
set of galax;ies in more than one filter. This subset could be 
obtained externally, from other existing or dedicated obser- 
vations by, for example, HST, or it could be obtained from 
within the Euclid-like survey itself, which we explore further 
here. We propose that with real data we would follow the 
procedure shown in Fig. [9] in which the parameters of each 
galaxy in the calibration sample are obtained from noisy im- 
ages of the galaxy observed in two filters. We assume that 
the galaxy model is known and use the estimated param- 
eters to simulate the galaxy as it would be observed in a 
single wide filter. Each galaxy is simulated in a ring-test to 
obtain the shear calibration bias. 
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Figure 8. Error on the mean multiplicative bias mi-p-^ ("^2;Fi 
omitted for clarity) as a function of number of galaxies in the 
calibration sample for a perfect calibration sample and assuming 
the distribution in mi-^p-^ obtained using the catalogue is Gaus- 
sian (blue crosses solid) and sampling from the true distribution 
(red circles). Results are also shown for a calibration sample im- 
aged in two- filters {Fi and F2) with a signal-to noise ratio of 100 
(green squares dashed), 50 (cyan stars dot-dashed) and 25 (ma- 
genta triangles dotted) in the _Fi filter. Approximate lower limit 
on the error on the mean bias for calibration using two filter imag- 
ing assuming a CWW-Sbc (black dashed) and CWW-Im (black 
dot-dashed) spectrum for the disk. Grey panels as in Fig. \3\ 
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Figure 9. Schematic of our proposal for how to find the 1 filter 
bias from calibration data observed in 2 filters. 



We first estimate the number of galaxies required to 
calibrate the bias from single broad-filter imaging assuming 
perfect information about the galaxy properties and spectra 
in the calibration sample (Section 16. ip . We then consider in 
Section 16.21 how well we can estimate the bias from single 
broad-filter imaging using information from noise-free im- 
ages of galaxies observed in the Fi and F2 filters. Finally, in 
Section [6.31 we investigate the effects of adding noise to the 
two-filter images of the calibration data. 

6.1 Idealised calibration data 

Here we estimate the number of galaxies required to cal- 
ibrate the bias on cosmic shear measurements from single 
broad-filter imaging. For a tomographic weak lensing analy- 
sis the galaxies will be divided up into ten or more redshift 
bins. To correct the bias we need only the mean bias on the 
shear at each redshift, and not the bias of each individual 
galaxy (see Appendix). 

Consider the optimistic case in which we have perfect 
information for each galaxy in a calibration sample, for ex- 
ample the spectrum at each point on each galaxy image. 
This information could be used to estimate the 1 filter shear 
measurement bias for each galaxy in the calibration sample. 
The average of these numbers could then be used to correct 
the shear power spectra measured from 1 filter imaging of a 



larger sample. Even though the information on each galaxy 
is perfect, this could only work if the calibration sample 
were a fair sample of the Universe, and contained enough 
galaxies. 

We use the distribution of bias values for single broad- 
filter imaging computed in Section [2 . 31 from the catalogue to 
represent the full population of galaxies in a given redshift 
bin. We draw rig random samples from the distribution in 
m to estimate the one sigma uncertainty on the mean mul- 
tiplicative shear bias for a sample of rig galaxies. The results 
are shown as a function of rig in Fig. |8] as red circles. 

We also show the results assuming the distribution in m 



is Gaussian, so that a{{m)) = a{m)/{ng 



1 



(blue line). 



The results are similar, showing that the distribution in m is 
well represented by a Gaussian. Therefore for a calibration 
sample in which the galaxy shapes and spectra are known 
perfectly we would need 3 x 10'^ galaxies to calibrate the 
bias on the shear to a factor of 10 better than the Euclid- 
like survey requirement. 

6.2 Two-filter calibration data 

We now consider how well the bias from single broad-filter 
imaging can be calibrated when the galaxies in the calibra- 
tion sample are imaged in just two filters. In this case we 
will not know the bulge and disk spectra precisely, but only 
linear fits to the bulge and disk fiuxes in each filter. There 
may also be some error on the galaxy shape parameters due 
to the approximate spectra. 

The procedure is outlined in the fiow diagram in Fig. 1 101 
We first simulate the 'average' galaxy from the catalogue 
(see Section 12. 4|) using a CWW-E spectrum for the bulge 
and either a CWW-Sbc or CWW-Im spectrum for the disk. 
We convolve the bulge with the bulge PSF and the disk with 
the disk PSF (i.e. simulate the true PSF-convolved galaxy 
image). The shape parameters of the galaxy and linear ap- 
proximations to the bulge and disk spectra are computed 
using two filter imaging in Fi and F2 as in Section [S] This 
is the information that we would obtain for a given galaxy 
in our calibration sample that could be used to estimate the 
bias on the shear for single broad-filter imaging. We thus 
simulate the galaxy using this information and perform a 
ring-test as in Section [4. II to obtain an estimate for the bias 
on the shear when the galaxy is imaged in the Fi filter only. 

We investigate the difference between the true mpi and 
rriFi measured from the calibration sample for a range of 
bulge-to-disk half-light radii and galaxy redshifts. The re- 
sults are shown in the right-hand panels of Fig. [7] for noise- 
free galaxy images. The error on from calibrating in 
two filters depends on the ratio of the bulge-to-disk half- 
light radii, the disk spectrum and the galaxy redshift. 

We estimate the mean calibration error (i.e. average dif- 
ference between the true and estimated mpi) over a popu- 
lation of galaxies. We create a population of galaxies at a 
given redshift with a CWW-Sbc spectrum for the disk, and 
the usual CWW-E spectrum for the bulge, and use the ratio 
of bulge to disk sizes from the catalogue. We calculate the 
error on mpi, as shown in the right panels of Fig. [T] For 
practical purposes we restrict our calculation to the range 
of values on Fig. [7] and if ro,b/?'c,d is less than the minimum 
(maximum) value on the a;-axis of the plot then the error is 
set equal to the value at the minimum (maximum) i-value 
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Figure 10. Flowchart to explain how we test whether two filters are sufficient to do the calibration. One filter shear calibration biases 
from the estimated galaxy parameters are compared to those from the true galaxy parameters. The estimated galaxy parameters include 
estimated spectra. 



shown. We repeat the calculation for each of the four redshift 
values plotted and estimate the mean error on mp^ over red- 
shift by weighting the contribution from galaxies at each red- 
shift u sing the galaxy redshift distribution from lSmail et al] 
||1994D : P{z) oc z°'exp[-(lAlz/ z^f], where a = 2, /3 = 1.5 
and the median redshift for a Euclid-like survey, ~ 0.9. 
We repeat the exercise using the CWW-Im spectrum for the 
disk. The error on niPi using the CWW-Sbc and CWW-Im 
spectra are shown by the lower dashed and upper dot-dashed 
horizontal hues respectively in Fig. [S] The error on the mea- 
sured bias cannot be reduced below the values shown by 
these horizontal lines by increasing the number of galaxies 
in the calibration sample. 



6.3 Noisy two-filter calibration data 

We repeat the calculations of the previous section but adding 
noise to the two-filter simulated images. The noise causes a 
spread of estimated galaxy parameters and thus increases 
the width of the distribution of m values and could induce 
an extra bias. 

A broader distribution of m values means that more 
galaxies will be needed in the calibration sample to measure 
the mean to the required accuracy. We add the variance 
of the m distribution obtained from the catalogue with the 
variance of the m distribution for the 'average' galaxy from 
the catalogue whose parameters were obtained from two- 
filter imaging of noisy images of the galaxy. The number 
of galaxies required in the calibration sample for different 
signal-to-noise values is shown in Fig. [8] The signal-to-noise 
is the value measured in the Fi filter. We assume the same 
exposure time in the Fi and F2 filters and so the signal-to- 
noise in F2 is less than the values quoted. Note for illustra- 
tion we assume that for a given galaxy the broadening of 
rnpi due to noise is independent of the galaxy shape param- 
eters and colour gradient and equal to the broadening for 
the 'average' galaxy. 

We check that the mean bias for single filter imaging 
in Fi estimated from two filter imaging of noisy calibration 
data is not wrong by more than the requirement on the 
bias. We find that the difference between the true bias from 
imaging in a single wide filter and the mean estimated bias 
over many noise realisations is at least a factor of two less 



than the Euclid requirement down to a signal-to-noise of 25. 
We note that there is no indication that the bias difference 
increases as the SNR decreases from 100 to 25. 



7 DISCUSSION 

To capitalise on gravitational lensing as a cosmological probe 
it can be attractive to use wide-filter imaging to gather as 
much light as possible and constrain cosmological parame- 
ters, including the dark energy equation of state parameters, 
to high precision. Such wide-filter experiments may however 
cause an unacceptable loss in accuracy. In this paper we con- 
sider the shear measurement biases arising from single filter 
imaging due to the wavelength-dependence of the PSF and 
the existence of colour gradients across galaxy images. 

The PSF is modelled as a Gaussian at each wave- 
length with ellipticity 0.05. We have assumed a particular 
wavelength dependence of the PSF size. In practice differ- 
ent telescope designs will have different dependencies de- 
pending mainly on the mirror size, the amount of jitter 
and the amoun t of c harge diffusion. As can be seen from 
ICvpriano et al. the optics dominate the PSF wave- 

length dependence and the biggest likely change would be 
the relative contribution of the other two (relatively wave- 
length independent) terms. A decrease in the jitter and/or 
amount of charge diffusion would increase the wavelength 
dependence of the PSF. In the worst possible case of zero 
charge diffusion and no jitter, the slope of the PSF size as 
a function of wavelength is approximately doubled, which 
we would naively expect to double the shear measurement 
biases. In practice the change would be smaller, and if the 
diffusion or jitter were larger than assumed then the biases 
would be smaller. We have assumed that the PSF elliptic- 
ity is independent of wavelength, but in practice additional 
components such as beam splitters could introduce wave- 
length dependent ellipticity. We assume that the ellipticity, 
orientation and wavelength-dependence of the size of the 
PSF are all known. 

We consider simple bulge plus disk galaxies modelled 
using two co-centered, co-elliptical Sersic profiles, each with 
a different spectrum. The bulge is modelled as a de Vau- 
couleurs profile and the disk as an exponential. We assume 
the galaxy model is known and fit for the centre, elliptic- 
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ity and orientation of the galaxy, the sizes of the bulge and 
disk components and the bulge-to-total flux ratio. We also 
assume the composite spectrum emitted by the galaxy is 
known from ground-based observations and do not take into 
account errors on the shape of the spectrum or total flux 
emitted. In practice the different passbands and apertures 
used from ground and space will make it difficult to ob- 
tain the true total flux which will make the estimated biases 
noisier and potentially biased. 

The results will depend on the complexity of the 
galaxy model. For example it is possible that non-elliptical 
isophotes of galaxies can produce non- negligible additive 
shear calibration biases (Hirata & Bernstein priv. comm.). 
However an investigation using more complex models than 
the elliptical isophote model considered here is beyond the 
scope of this work. In particular we note that extending the 
calculation to more complex galaxy models is non-trivial 
because current shape measurement methods do not meet 
the accuracy requirements for a Euclid-like survey for real- 
istic signal -to-noise galaxies with radially- varying ellipticity 
isophotes (|Bridle et al.ll201ol ): though considerable progress 
has been made on un derstanding the limit ations of ex - 
isting methods l|Voigt fc Bridle _201P; Mclchior et al.ll201C ') 
and on the developm ent of new methods l|Bernsteinl 201C ; 
iKitching et"al]|2010t ). 

We first measure the bias on the shear if we ignore 
colour gradients and image galaxies in a single filter. In the 
simulations of the 'true' galaxy images we model the bulge 
spectrum using a CWW-E spectrum and the disk using ei- 
ther a CWW-Sbc or CWW-Im spectrum. We do not take 
into account evolution of these spectra with redshift as this 
is much smaller than the difference between the CWW-Sbc 
and CWW-Im spectra and we show results for both. In the 
fits to the simulated images we assume the spectrum at each 
point in the galaxy image is equal to the composite spectrum 
from the bulge plus disk emission. 

We measure the dependence of the bias on the width of 
the filter for example galaxy parameters. The bias drops by 
a factor of approximately 5 when the filter width is halved. 
The bias depends strongly on the galaxy shape parameters 
as well as on the colour gradient. For a 'typical' galaxy and a 
filter width of 350 nm the bias is a factor of 3 below the Eu- 
clid-like survey systematic equals random error requirement 
of 10~^, but the bias can be more than a factor of 10 above 
the requirement for some galaxy properties considered. 

We show that, to first order, tomographic cosmic shear 
two point statistics depend on the mean shear bias over 
the galaxy population at a given redshift. We compute the 
mean bias for single broad-filter imaging, averaging over a 
catalogue of galaxies from a published study of HST ob- 
servations of galax ies in the Keck LRIS DEEPl survey 
l|Simard et al.ll200^ V The catalogue is used to represent the 
range of galaxy colour gradients, shape parameter ratios 
7'e,b/7'e,d and bulge-to-total flux ratios in the Universe. We 
caution, however, that the catalogue has a complicated se- 
lection function (see Section 12. 3p . and that this must be 
taken into account when drawing conclusions for future cos- 
mic shear surveys. We attempt to investigate the sensitivity 
of the results to the catalogue by plotting the bias depen- 
dence on different galaxy parameters, and by re-computing 
the mean bias for different colour cuts on the catalogue. 

We find that the mean bias over the full catalogue is 



approximately a factor of 3 below the Euclid-like survey re- 
quirement. (The mean bias over the catalogue is similar to 
the bias for the fiducial galaxy because the parameters of the 
fiducial galaxy are chosen to represent the 'average' galaxy 
from the catalogue). However, to ensure the total bias arising 
from all shear measurement systematics is below the Euclid- 
like requirement, the mean bias from the colour dependence 
of the PSF should be of order a factor of 10 below the sys- 
tematic equals random error requirement. We thus conclude 
in Section |4] that single broad-filter imaging in a future Eu- 
clid-like survey may not allow the accuracy requirements on 
the bias to be met, depending on how accurately the cata- 
logue from lSimard et all (|2002r ) represents the cosmic shear 
survey being considered. 

We assume that the mean bias obtained using the full 
catalogue is optimistic because we find it may be up to a 
factor of 2 larger than 10"'^ if we modify the galaxy red- 
shift distribution by making colour cuts on the catalogue. 
However, it is possible that the mean bias computed using 
galaxy parameter distributions from the catalogue is larger 
than the true mean bias in a Euclid-like survey. This is be- 
cause the photomet ric and structural parameters obtained 
in the lSimard et al.l (j2002) study will not precisely represent 
the underlying distribution of colours and shapes of galax- 
ies in the DEEPl survey a s a result of fitting a simplified 
model to noisy galaxies fsee lHaussler et al ] |2007D . It is likely 
that this will result in extreme bulge-to-total fiux ratios and 
colours, increasing the mean bias. 

The bias from single broad-filter imaging can be reduced 
to a negligible level by narrowing the filter, however the sur- 
vey time or efficiency would have to increase to maintain 
the same cosmological constraining power. We also consider 
two other possible options: (i) multiple-filter imaging of the 
whole survey and (ii) calibrating the bias from single broad- 
filter imaging on a sub-set of galaxy images observed in more 
than one filter. We first investigate how much the biases may 
be reduced by imaging in two filters. We consider an addi- 
tional narrower filter with width equal to half the width of 
the wide filter and fully contained within the wide filter. The 
colour information available from a two-filter configuration 
allows us to determine linear approximations to both the 
bulge and disk spectra across the wide filter. Using the bulge 
and disk PSF models calculated from these linear spectra, 
and including the information from the composite spectrum, 
the biases are reduced by an order of magnitude below the 
biases from single broad-filter imaging. Thus, for the simple 
bulge plus disk galaxy model considered, only one narrow 
filter is needed in addition to the wide filter to remove the 
bias from the PSF wavelength dependence to the required 
accuracy. We note that if the galaxy has more than two 
components, each with a different spectrum, then two-filter 
imaging is sufficient to estimate linear approximations to the 
spectra of each of the components, provided the components 
are resolved. 

We also consider the possibility of calibrating the bias 
from single broad-filter imaging on a sub-set of galaxy im- 
ages. For a tomographic cosmic shear analysis we need the 
mean multiplicative bias for each redshift bin to be mea- 
sured to an accuracy of 10~^ (i.e. a factor of ten below the 
Euclid requirement). We note that we expect the mean bias 
to be redshift-dependent because galaxy morphologies and 
spectra evolve with redshift. We first compute the number of 
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galaxies needed in the calibration sample to calibrate a Eu- 
clid-like survey assuming perfect knowledge of the galaxy 
shape parameters and spectra. Using the distribution of 
galaxy shapes and colour gradients in the catalogue as a 
representation of the population of galaxies in the Universe 
we find that we need approximately 3 x 10"^ galaxies per 
redshift bin. 

We then investigate whether the information available 
from noisy two-filter imaging is sufficient to calibrate the 
bias from single broad-filter imaging. We assume a constant 
exposure time in the two filters so that the signal-to-noise 
in the narrower filter is lower than in the wide filter. We 
find that for a signal-to-noise in the wide filter equal to 25 
the number of galaxies required to calibrate the bias in each 
tomographic redshift bin is of order fO'*. 

We caution, however, that the calibration strategy as- 
sumes the galaxy model is known and further investigation 
is required to determine whether the procedure will work 
on low signal-to-noise images of real galaxies with complex 
morphologies. We check that the error on the calibration 
bias measured from two-filter imaging of noisy galaxies is 
sufficiently small down to a signal-to-noise of 25. We do not 
check lower signal-to-noise values due to computational lim- 
itations. We note that in practice it may be possible to stack 
images with similar properties, e.g. colours, within each red- 
shift bin to reduce the noise. 

We have used a simple galaxy model in our proposed 
calibration strategy because this may capture the main 
colour gradient effect that causes the shear measurement 
bias. Ifowever this may be insufficiently flexible in practice, 
depending on the complexity of galaxies in the Universe and 
any additional shear measurement bias this causes. For this 
simple galaxy model we find that the shear measurement 
bias is relatively insensitive to the true galaxy ellipticity 
(lower left panel of Fig U), and therefore galaxies from a 
calibration sample could be imaged with lower PSF qual- 
ity than for the main cosmic shear sample. Qualitatively 
this is because the calibration sample provides information 
on galaxy colour gradients which has a less stringent image 
quality requirement than the very stringent requirement for 
cosmic shear. 

We note that the bias value we obtain from the galaxy 
catalogue cannot be taken to be the true bias for a Eu- 
clid-like survey because it contains an insufficient number of 
galaxies and the selection criteria are not matched. The sam- 
ple contains 632 galaxies over a range of redshifts whereas 
we conclude that a calibration sample of around f 0* galaxies 
per tomographic redshift bin is required. Our conclusions on 
the number of galaxies needed in a calibration sample should 
be more robust though, because the catalogue gives us in- 
formation about the range of galaxy types in the Universe. 

We note that around rb = rd the multiplicative bias 
changes sign. We expect the bias to be zero at rb = rd 
for the simple case where the spectral index of the bulge 
and disk are the same, because the radial intensity profile 
would be identical for each component and thus the true 
PSF-convolved galaxy image equal to the true galaxy image 
convolved with the PSF for the composite spectrum. We find 
that when the bulge and disk spectral indices are different 
the bias is negligible for rb/rd close to f. We note that the 
average of the Simard catalogue occurs near to this point 
(rb/rd = l.f) and therefore this contributes significantly to 



the small bias obtained due to the cancellation of large biases 
of either sign. Again, if the Simard catalogue is atypical then 
this could lead to a larger bias from realistic data. 
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APPENDIX: IMPACT ON SHEAR POWER 
SPECTRA OF DENSITY DEPENDENT SHEAR 
CALIBRATION 

Galaxies with a particular colour gradient will tend to be 
clustered, for example clusters of galaxies contain a lot of 
ellipticals which tend to have a small colour gradient, and 
the converse will be true in voids. In this Appendix wc check 
how this will affect the tomographic lensing power syjcctr um. 

This calculation is related to that bv iCuzik fc BernsteirJ 

l|20Q5l ) who calculate the effect of spatially varying mul- 
tiplicative shear calibration biases which are uncorrelated 
with the shear. Here we assume that the multiplica- 
tive bias is correlated with the density perturbations and 
therefore to the shear. There are some parallels between 
this Appendix and calculations of the effect of galaxy 
number density wei ghted shear measuremen ts (e.g. see 



Hhata & Seliak 2004) or s ource-lens clustering ([Bernardeaul 
1998 ; Hamana et al. 20o2; iForero-Romero et allBoOTi). and 



the s hear-intrinsic alignment or GI effect ( Hirata fc SeliakI 
l2004h . However, the former calculations do not consider 
shear calibration bias. 

In this derivation we assume that galaxies have a 
stochastic multiplicative shear measurement bias m so that 
the observed shear d is related to the true shear 7 by 



d = (1 -I- m)7 



(18) 



where m is drawn from a probability distribution whose 
shape depends on the spatial location Pi{m\0, z). Here we 
assume that the whole probability distribution is stretched 
or squeezed uniformly according to the local matter density 



Pr(mj<5) = f{{l + bmS)m) 



(19) 



where bm quantifies the sensitivity of the multiplicative bias 
to the mass distribution. We would expect bm to be neg- 
ative to obtain the expected reduced multiplicative bias in 
a local overdensity. Here we choose to assume bm is inde- 
pendent of scale but is an arbitrary function of redshift 
b,n{z), which would depend on galaxy evolution. It would 
be possible to generalise this derivation for density indepen- 
dent s hear calibration biases discussed in lGuzik fc Bernstein! 
(|2005l ) by adding in terms to the above equations, but we 
choose to focus the discussion here on the new effect. 
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For first order calculations of the shear two point statis- 
tics the relevant property of this distribution will be the 
mean, fa, which can be written 



fh{0, 2) = (1 + bra{z)S{e, z))fh{z) 



(20) 



where m is the mean multiplicative shear calibration bias 
averaged over the sky at a given redshift, since (5) = 0. The 
description of the probability distribution could therefore be 
more general than described in the paragraph above. Note 
that to achieve the usual effect of decreased biases at high 
densities, hm needs to be negative. At very high densities this 
equation predicts very negative multiplicative shear biases, 
which would fit with very blue cores in galaxies in the centres 
of clusters. This is not particularly physical and therefore 
the model should be improved by using a more sophisticated 
function of 5 which could affect the cancellations in the final 
result. For example some inverse function of (1 + 5) could 
have a more realistic effect at high densities. 

The main cosmic shear statistic of interest for future 
surveys is the tomographic shear power spectrum, which 
may be defined in terms of the true shear 7 

C«f = {{liillj)4,i)i 

where the inner angle bracket denotes an average over direc- 
tion of the vector I and the outer angle bracket over reali- 
sations. We use G to denote true gravitational shear. Here i 
and j denote the index of the tomographic redshift bins used 
for the correlation. Tilde denotes the usual Fourier transform 
in the plane of the sky 



(22) 



The shear 7 is related to the mass density 5 by the usual 
equations 



7(0,2) = D{e) * K{e,z) 

where 

k{0,z)= / dz' q{z, z')S{0, z') 



(23) 
(24) 



using the Limber approximation and writing the three di- 
mensional mass density per unit comoving volume in terms 
of position on the sky and redshift. q is the usual lens weight 
function but the only relevant property of q for this paper 
is that q{z, z') is zero when z — z . 

We assume that tomographic redshift bins are prede- 
termined according to some observable property which is 
independent of the density field so that the number density 
of galaxies in redshift bin i as a function of position on the 
sky and redshift ni(Q, z) may be separated out as 



n.(6». 



n,(z)(l + (55(6»,2)) 



(25) 



where ruiz') is the mean number density of galaxies at red- 
shift z in redshift bin i averaged over a patch of sky large 
enough that the mean galaxy density fluctuation 5g(0, 2) is 
zero, normalised such that 



{z)dz = 1. 



(26) 



We can only calculate statistics of the observed shear d 
and we here assume this is done by producing a map of d 
on the sky for each redshift bin. We assume this is done by 



averaging observed shears in pixels on the sky, therefore this 
map can be written in the limit of infinite source density as 

d^(Q)= dz dm n,{0,z)Pr{m\e,z)d{m,e,z) (27) 



which can be written to first order in the small quantities as 
(6>) = 7, (6>) 79, (6») + 7m, (6>) (28) 
where we have defined 



7,(6>) = dz n,{z)-y{0,z) 



dz ni{z)Sg{0, 2)7(0, 2) 



(29) 



(30) 



7m(0) = / dz ni{z) / dmPr{m\0, z)mdm~f{0, z) (31) 



(32) 



dz ni{z)m{0, 2)7(0, 2). 



We will measure the power spectrum of the observed 
shear d 



(21) = {{dudij)ci,e)r 



(33) 



which can be written to first order in the perturbative quan- 
tities as 



where we have defined 



C\ 



Gg 
lij 

Gm 
lij 



(34) 

(35) 
(36) 



where is the same as the term written with the same 
notation in Ijoachimi fc Bridle! l|2010l ) but C«™ is different 
because in this paper we use m to denote the multiplicative 
shear calibration bias whereas in that work m was used to 
denote magnification which we ignore here. 

The first term above is the usual shear term and 
the only one usually considered, the second takes into ac- 
count the fact that we measure the galaxy number density 
weighted shear (which contributes to source-lens clustering) 
and the third term takes into account the possible depen- 
dence of the measured shear on a multiplicative shear mea- 
surement calibration bias. We write down the source-lens 
clustering term here because the effect may at first sight 
seem correlated with the effect discussed in this paper, how- 
ever from the above equations it can be seen that to the level 
of approximation considered here the terms are independent, 
and we thus discuss it no further. Note that the final two 
terms can be non-zero for overlapping redshift bins. 

The cross shear-multiplicative bias term may be further 
expanded as 



C, 



Gmb 



where we have defined for convenience 



1fhi{0) 



7mM(^) = 



Gm, 
■■lij 

CGmh 
lij 



dz ni{z)m{z)j{0, z) 

dz ni{z)rh{z)b,niz)S{0, z)'y{0, z) 

{{lul'kij)<t>i)i 

{{lli'imblj)<t>l)r- 



(37) 



(38) 

(39) 

(40) 
(41) 
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We now focus on the second term above and substitute 
in the relation between shear and mass density. We find 

'Sihii = J dz ni{z)fh{z)hm{z) j dz' q{z' , z)Ae{z, z') (42) 

where we have defined 

Ae{z,z') = Se{z)*{DA{^')). (43) 
Therefore 

C«™'' = J dzgi(z) j dz nj(z')ih{z')h,n{z') j dz" q(z" , z')Bd,i{z, z , z"){4:4) 
where we defined 

Bd*i{z,z',z") = {{DMz) {W) * {Di5i{z)))*)^i)r. (45) 

If we assume that mass at one redshift is uncorrelated 
with mass at another redshift then Bo*e{z, z' , z")) is zero 
unless z = z' = z" . Taken together with the fact that 
q{z,z') = when z = z' it would follow that Cf^^ = 
at all redshifts. In any case this term depends on the three 
point function of the mass distribution and is expected to 
be even smaller than the source-lens clustering term due to 
the additional factor of m. 

The remaining term C^™ is much simpler and can be 
written 



Cf^r = j dzg,{z)g^,{z)P{k = £/x, z) (46) 
where we have defined 

gmi{z)= j dz'ni{z)rh{z')q{z,z'), (47) 



X is the comoving distance to redshift z and k is the comov- 
ing wavenumber. In the special case where m is independent 
of redshift then 

= fhC^f (48) 

as expected for the simple case where d = (1 + fh)"f (note 
the additional term CYlj^ which makes up the factor of 2). 

Therefore we have found that the cosmic shear power 
spectrum is modified to first order by a term proportional 
to the average of the multiplicative shear calibration bias. 
In detail there is an additional power spectrum which is 
identical to the usual lensing power spectrum but in which 
the lensing weight function is weighted by the mean multi- 
plicative shear calibration bias m at each redshift, spatially 
averaged over the sky. This is why we require for calibra- 
tion of this effect the mean multiplicative shear calibration 
factor at each redshift. This should be spatially averaged 
rather than galaxy number density averaged. The number 
of redshift bins in which m should be calculated depends on 
how smoothly it changes with redshift. 

(49) 
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